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A set of exact quasi-local conservation equations is obtained in the (1+1) -dimensional descrip- 
^v^j ", tion of the Einstein's equations of (3+l)-dimensional spacetimes. These equations are interpreted 

■ as quasi-local energy, linear momentum, and angular momentum conservation equations. In the 

asymptotic region of asymptotically flat spacetimes, it is shown that these quasi-local conservation 
£\j , equations reduce to the conservation equations of Bondi energy, linear momentum, and angular 

momentum, respectively. When restricted to the quasi-local horizon of a generic spacetime, which 
is defined without referring to the infinity, the quasi-local conservation equations coincide with the 

Q conservation equations on the stretched horizon studied by Price and Thorne. All of these quasi- 
local quantities are expressed as invariant two-surface integrals, and geometrical interpretations in 
• terms of the area of a given two-surface and a pair of null vector fields orthogonal to that surface 

are given. 

(N ' I. INTRODUCTION AND KINEMATICS 

-xf : 

, For the past few decades, there has been enormous progress in general relativity since the pioneering works of 
the late Professor A. Lichnerowicz on mathematical relativity. His contributions to general relativity are diverse as 
£S) \ well as profound, not least because he put the Einstein's equations on a firm mathematical foundation as partial 
differential equations and theories of connections 0. Connections also plays important roles in Yang-Mills gauge 
I theories, since gauge theories are nothing but theories of connections coupled to matter fields. Therefore, in this 
£jl International Conference commemorating Professor A. Lichnerowicz, it seems appropriate to discuss a relatively 
^ 1 unknown formalism of general relativity, which is based on the very idea of connections. 

(3JT), This note is about (l+l)-dimensional description of general relativity of (3+l)-dimensional spacetimes, treating 
the remaining 2-dimensional spatial dimensions as a fibre space. In this framework, all the notions in the theory 
. — of fibre bundles such as a fibre space, connections, and the structure group appear naturally. Instead of going into 
the details of the formalism itself^], |^, ||, however, I will describe the key ideas briefly, mainly to fix the notations, 
and then quickly move on to discuss issues that arc more immediate, namely, the problem of defining the quasi-local 
conservation equations using the (l+l)-dimcnsional formalism of (3+l)-dimcnsional spacetimes. 

Let us begin by mentioning a few facts about quasi- local conservation equations. In general relativity there have 
been many attempts to obtain quasi-local conservation equations^, p|, |6|, f^, p[ p| fLof. One of the motivations of these 
efforts is the expectation that quasi-local conservation equations allow us to predict certain aspects of a quasi-local 
region of a given spacetime without actually solving the Einstein's equations for that region. Recall that in the 
Newtonian theory, the conservation of total momentum immediately follows from Newton's third law, 

Ftotal = ^£>) = 0, (1) 

i 

which is no more than the consistency condition implementing Newton's second law. In general relativity, the consis- 
tency conditions for evolution are already incorporated into the Einstein's equations through the constraint equations, 
from which global conservation equations were found. In this note we will show that, from the Einstein's equations 
in the (l+l)-dimensional description, one can find conservation equations of a stronger form, namely, quasi-local con- 
servation equations. These equations, which are integro-differential equations over a compact two-dimensional space, 
are naturally interpreted as quasi-local energy, linear momentum, and angular momentum conservation equations fli"l|. 
Let us consider the following line element 

als 2 = -2dualv - 2hdu 2 + c a p ab (dy a + A^du + AJ*dv) (dy b + A*du + A*dv) , (2) 

where +, — stands for u, v, respectively || [|, |TJ, |TJ, |TJ, |TJ, [TtJ . To understand the geometry of this metric, it is 
convenient to introduce the following vector fields, 

8+ := d+ - A -8 a , (3) 
d- := cL - A*d a , (4) 
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where we defined the following short-hand notations 

d „ d 



d+:= d^> d - := d^r ^ : -^( fl - 2 ' 3 )- ( 5 ) 

The inner products of the vector fields {<9±, d a } are given by 

< d+, d+ >= -2h, <<§+,<§_>= -1, <d-,d->=0, 

< d±, d a >= 0, < d a , d b >= e a Pah . (6) 

The hypersurface u — constant is a null hypersurface generated by the out-going null vector field <9_, which is 
orthogonal to the vector fields {d a }- Notice that v is the affine parameter of the out-going null vector field. The 
hypersurface v = constant is generated by the vector field d+ whose norm is —2h, which can be either negative, zero, 
or positive. The intersection of two hypersurfaces u,v = constant defines a spacelike compact two-surface N2, which 
are coordinatized by y a . The metric on N% is decomposed into the area element e CT and the conformal two-metric p ab , 
which is normalized to have a unit determinant 

dct p ab = 1. (7) 

In the terminology of the fibre bundles, the base manifold is the (l+l)-dimensional space-time coordinatized by (u, v), 
and the fibre space is 2-dimension spacelike space N2. The vector fields {d±}, which are orthogonal to {d a }, is the 
horizontal vector field, and {d a } is tangent to the fibre space A^. The fields A± are the corresponding connections 
valued in the diffeomorphisms of the two-surface N2 |^, 0] . 

For later uses, we shall write down the future-directed in-going null vector field n and out-going null vector field /, 
orthogonal to two-surface A^ 2 at each spacetime point. They are given by 

n:=d+-hd-, (8) 

/:=£_, (9) 

and are normalized such that 

< n, I >= -1. (10) 

If we further assume that A_ a = 0, then the metric (||) becomes identical to the metric studied in fl6| . In this note, 
however, we shall retain the field, since its presence will make the A^-diffeomorphism invariant Yang-Mills type 
gauge theory aspect of this formalism transparent. Apart from the iV2-diffcomorphism invariance, there are other 
residual symmetries that preserve the metric (^|), which are the reparametrization of u, and the transformation that 
shifts of the origin of the affine parameter v at each point of N2 [0 . 

The complete set of the vacuum Einstein's equations are found to be|lq| 



(a) e a D + D-a + c a D-D + a + 2c CT \D + <t)(D-<t) - 2e ff (D_ h) (L>_ a) - -c 2a p ab F + ^F + b 



+c a R 2 - /ie CT {(D_a) 2 - ^p ab p cd (D_p ac )(D_p bd )} = 0, (11) 

(b) -e a D\a - ^e a (D+a) 2 - e° (D_h)(D+a) + e a (D + h)(D-a) + 2hc a {D^h){D-<j) 
+e°F + ^d a h - -^p ab p cd {D +Pac ){D + p bd ) + d a (p ab d b h) 

+h[ - <F{D+a){D-a) + ^ p ab p cd (D + p ac )(D_p bd ) + ^e 2 ° Pab F + lF + l - c a R 2 ) 
+h 2 e°[{D-af - ±p ab p cd (D_p ac )(D_ Pbd )} = 0, (12) 

(c) 2e°(D 2 _a) + c ^ D _ a f + ^p ab p cd (D^p ac )(D^p bd ) = 0, (13) 

(d) D-(e 2 °p ab F + b ) -e°d a {D_o) - ^ p bc p de (D_p bd )(d aPce ) + d b (e° p bc D_p ac ) = 0, (14) 

(e) -D + (e 2 ° Pab F + b ) - c°d a {D+a) - ^ p bc p de (D +Pbd )(d a p ce ) + 8 b (e° p bc D +Pac ) 
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+2he a d a {D_a) + he° p bc p de (D_p bd )(d a p ce ) + 2e°d a (D_h) - 2d b (he° p bc D_p ac )j = 0, (15) 
(/) -2e a D 2 _h - 2e CT (D_/i) (£>_+) + e a D + D-<r + c a D-D + a + e° (D + a)(D_a) 
+ \z a p ab p cd (D + p ac ){D_p bd ) + c 2a Pab F + «F + b - 2he°{D 2 _a + \{D-a) 2 

+ \p ab p cd {D„p ac ){D_p bd ))=Q, (16) 
(g) h[e a D 2 _ Pab - e a p cd (D-p ac )(D-p bd ) + e° (D- Pab )(D-cr)} 

~e° (D + D_p ab + £>_D+p a6 ) + ^e a p cd {{D- Pac )(D +Pbd ) + (D_p bc ){D + p ad )} 
~e a {(D_p ab )(D+a) + (D +Pab )(D_<j)} 

+e a {D_ Pab )(D_h) + ^e 2 °p acPbd F + lF + d - -U 2a PabPcd F + c lF + d = 0. (17) 

Here R 2 is the scalar curvature of iVjj, and we defined the diff^-covariant derivatives as follows, 

F + « := d+A_ a - d-A« - L4+, (18) 

D±a:=d±<T-[A±,<T] L , (19) 

D±h:=d±h-[A±,h] h , (20) 

D±Pab ■= d± Pab - [A±, p\ Lab . (21) 

The bracket [A±, f]h a b--- is the Lie derivative of f ab - - along the vector field A± := A±d a , defined as 

[A±, f] Lab ... := A£d c f ab ... + f cb ...O a A^ + hc-dbAl - w(d c A^)f ab ..., (22) 

where w is the weight of the tensor density f ab .... One can also compute the scalar curvature R of the metric (^) and 
integrate it over spacetime. It is given by 



la = Jdudvd 2 ye cr R 

= Jdudv d 2 y Lq + surface integrals, (23) 
where the "Lagrangian" function Lq is given by|| 



+- r +- 



L° = -^ 2a p ab F + «F + b + e°(D + a)(D_a) - -e" P ab P cd (D +Pac )(D_ Pbd ) - e° R 2 
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-2e a (D_h)(D_a) - he a {D_<y) 2 + -he* p ab p cd (D_p ac )(D_ Pbd ). (24) 

One can easily recognize that this "Lagrangian" function Lq is in a form of a (l+l)-dimensional field theory 
Lagrangian. In geometrical terms the function Lq describes how the (l+l)-dimcnsionsal spacetime and 2-dimensional 
fibre space are imbedded into an enveloping (3+l)-dimensional spacetime. Each term in (Q) is manifestly diffA^- 
invariant, and the y a -dependence of each term is completely "hidden" in the Lie derivatives. In this sense we may 
regard the fibre space N 2 as a kind of "internal" space as in a Yang-Mills theory, with the infinite dimensional group 
of diffeomorphism of N 2 as the Yang-Mills gauge symmetry. Thus, the above function Lq is describable as a (1+1)- 
dimensional Yang-Mills type gauge theory interacting with (l+l)-dimensional scalar fields and non- linear sigma fields 
of generic types. 



II. A SET OF QUASI-LOCAL CONSERVATION EQUATIONS 

Notice that the four equations (|lT|) , (]l^) and (JT^) are partial differential equations that are first- order in £)_ 
derivatives. Therefore it is of particular interest to study these four equations, since they are close analogues to 
the Einstein's constraint equations in the usual (3+1) formalism. Thus, in this formalism, the natural vector field 



that defines the evolution is D-. Then the momenta ttj — {TT} l ,ir r7 ,ir a ,ir ab } conjugate to the configuration variables 
q 1 = {h, a, A? , p^} are defined as 



They are found to be 



9L . . 

= WW)' (25) 



7T ft = -2c ff (£_a), (26) 

tt ct = -2c a (D-h) - 2he a (D-o) + c a (D + a), (27) 

7T a = G 2a Pab F + h _, (28) 

n ab = he°p ac p bd {D_p cd ) - ^p ac p bd (D +Pcd ). (29) 



Notice that ir ab is traceless 



n a a = 0, (30) 
due to the identities that are direct consequences of the condition ([?]) , 

p ab D ±Pab = 0. (31) 

The "Hamiltonian" function Hq defined as 

H := niD^q 1 - L (32) 



is found to be 



where H is given by 



Hq = H + total divergences, (33) 



H = —e-^TTh + \he-^ 2 h - i e - 2CT p ab 7r a 7r 6 + p acPbd -K ab Tr ed 

+^ h (D+a) + ^ir ab (D +Pab ) + ±-e° p ab p cd (D + p ac )(D +Pbd ) + e CT i? 2 . (34) 

In terms of these canonical variables {717, q 1 }, the first-order equations (fill), (|l2|), and ( |l5| ) can be written as, after a 
little algebra, 



(i) TT ab D +Pab + 7r a D+a - hD+7T h - d+ [h ir h + 2e a D+a ] 
+d a (h Tr h A_? + 2A^c a D+a + 2he~ a p ab ir b + 2p ab d b h) ) = 0, (35) 

(ii) H - d+n h + d a (a^-kh + e~ V^fc) = 0, (36) 

(iii) d+n a - d b (A b n a ) - n b d a A b - TT a d a (J + d a n a - n h d a h - ir bc d a p bc 

+d b (ir bc p ac ) + d c (ir bc Pab ) - d a (iT bc p bc ) = 0. (37) 

These are the four first-order equations in the gauge (||), and it is these equations that we are concerned with in this 
note. Notice that the equations (35) and ( |36|) are divergence-type equations. If we contract the equations (^) by an 
arbitrary function £ a of {v, y b } such that 

d+e = 0, (38) 
then the resulting equation is also a divergence-type equation, 

7T ab £ S Pab + 1T*£ S <T + 7T h £^h + TT a £ S A^ - 0+ (£°7r«) + d a ( - £<V<r + T p bc + A^^b) = 0, (39) 

where £^ is the Lie derivative along the vector field £ := £, a d a - 
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The integrals of these equations over a compact two-surface N2 become, after the normalization by I/I671", 

T^U{u, v) = ^ jd 2 y (n ab D +Pab + it a D+a - , (40) 

^-P(u,v) = ^-L 2 yH, (41) 



d_ 

du' 



-L{u,v ] i) = ^jd 2 y(n ab £ iPab + ii^ i a-h£^ h - A^f^ (d+f = 0), (42) 
where in the last integral we used the fact that 

>d 2 y£sf=(ld 2 yd a (ef)=0 (43) 



for a scalar density / with the weight —1. Here U{u,v), P(u,v), and L(it, u;£) are invariant two-surface integrals 
defined as 

U(u, v) := ^ j d 2 y (h ir h + 2c° D+a) + U, (44) 
P(u,v):=-L jd 2 y{-K h ) + P, (45) 
L(u,v;0^^fd 2 y(en a )+L (8+^ = 0), (46) 
where C/, P, and L are undetermined subtraction terms. Notice that these subraction terms must be u-independcnt, 

du du du 

in order to satisfy the equations (^0|), ([4l|), and (|42"|), respectively. In general the subtraction terms are not unique, 
and the "right" subtraction term may not even exist at all in a generic situation. One natural criterion for the "right" 
choice of subtraction term would be that it must be chosen such that the quasi-local physical quantities reproduce 
"standard" values in the well-known limiting cases. 

One can write the r.h.s. of the equation (|40|) in a more symmetric and suggestive form as follows. To do this, let 
us contract the equation ([37]) with and integrate over N2 to obtain the following equation 

d 2 y (A + Q a+7r a ) - jd 2 y (ir ab £ A+ Pab + ^ a £ A+ a - h£ A+ Tr h S ) . (48) 

If we use the definition of diff7V 2 -covariant derivatives D± and the equation (f48|), then the equation (|o|) can be written 
as 

■^-U(u, v) = ^ jd 2 y (Tr ab d +Pab + iTvd+a - hd + 7r h - A + a a +7 r a ) , (49) 

where the integrand on the r.h.s. assumes the canonical form of energy-flux, which is typically given by 

To+^Y^nd^, (50) 

i 

where <j) 1 is a generic field and iti is its conjugate momentum. Notice that the r.h.s. of the conservation equations 
( [42| ) and (^) match exactly, if we interchange the derivatives in the integrands 

£ s < — > d+. (51) 

In a region of a spacetime where d/du is timelike, these quasi-local equations becomes quasi-local conservation 
equations, which relate the instantaneous rates of changes of two-surface integrals at a given u-time to the associated 
net flux integrals. Let us remark that, unlike the Tamburino-Winicour's quasi-local conservation equations Jll| which 
are "weak" conservation equations since the Ricci flat conditions (i.e. the full vacuum Einstein's equations) were 
assumed in their derivation, our quasi-local conservation equations are "strong" conservation equations since only the 
four first-order equations were used in the derivation. 
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It is interesting to notice that we can obtain yet another quasi-local conservation equation. This is simply achieved 
by writing the equation (Esh as 



d_ 

du 



j>d 2 y {A«v a ) = j>d 2 y (n ab £ A+ p ab + n <T £ A+ a - h£ A+ K h + ir a d + A^ , (52) 



which relates the instantaneous u-derivative of the two-surface integral on the l.h.s. to the net flux integral on the 
right. However, the r.h.s. of this equation is not quite "canonical" due to the last term. If we restrict the field A+ 
such that it satisfies the u-independent condition 

d+A£ = 0, (53) 

which is essentially the same condition d38| ) that £ a satisfies, then the last term in the r.h.s. of the equation ( |5^ ) 
drops out, and we obtain the following equation 

-^J(u,v) = j)d 2 y (^T ab £A + Pab + ir a £A + <? - h£A + ^hj ■ (54) 

Here J(u, v) is defined as 

J(u,v)~j^Id 2 y(A°n a ) + J (d+A« = 0), (55) 

where J is an undetermined subtraction term. The r.h.s. of the equation ( [54] ) now represents a flux of the canonical 
form 

Y,n£A + 4>\ (56) 

i 

just as the r.h.s. of the equations (M2) and (E9h do. 



III. GEOMETRICAL INTERPRETATIONS 

Remarkably, the two-surface integrals (p4|), (|45|), (p6|), and (p5|), which were derived using the metric (||), can be 
expressed geometrically, in terms of the area of the two-surface and a pair of in-going and out-going null vector fields 
orthogonal to that surface. In order to show this, we need to invoke the definitions of in-going and out-going null 
vector fields {n, 1} defined in the section |[ 

A. Quasi-local energy 

Let us first observe that the integral in (|44|) can be written as the Lie derivative of the scalar density e CT along the 
in-going null vector field n, 

<j> d 2 y (hTT h + 2c a D + a^j = 2 <j> d 2 y e° '(d+ct - hD^a^j 

- 2<fd 2 y£ n e' J . (57) 



One finds that 

j> d 2 y £ n e' T = £ n A, (58) 

where A is the area of N2 , 

A = <fd 2 ye a . (59) 
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The identity (58) follows trivially from the observation that the null vector field n is out of (in fact, orthogonal 
to) the two-surface, which means that the order of the integration over d 2 y and the Lie derivative £ n in ( p)7| ) is 
interchangeable. Thus we have 



hw h + 2c a D+a) = ^-£„A. (60) 

07T 



For a reference term U, let us choose 



U:=~£ n A, (61) 

07T 



where n is a future-directed in-going null vector field of a background reference spacetime ds 2 into which the two- 
surface N 2 (with the same metric e a p a b) is embedded, 

ds 2 = -2dudv - 2hdu 2 + e a p ab (dy a + A_?du + A^dv) (dy b + A^du + A^dv) . (62) 

Notice that n, which is given by 

is a function of h and A± only, the embedding degrees of freedom of the two-surface. Thus, the quasi-local energy 
of a given two-surface N 2 is defined relative to some fixed background reference spacetime, and is zero when the 
two-surface under consideration is embedded into the fixed background reference spacetime, i.e. when 

n — n. (64) 

Therefore, the quantity U(u, v), which will be interpreted as the quasi-local energy, becomes 

U(u,v):=^-£ n A-^-£ n A. (65) 

It is given by the rate of change of the area of a given two-surface along the future-directed in-going null vector field 
n, relative to some background null vector field n. Notice that this definition is entirely geometrical, referring to the 
area of a given two-surface and the orthogonal null vector fields n (and n) only. 



B. Quasi- local linear momentum 

The two-surface integral ( J45| ) can be also written geometrically in a similar way. It becomes 

^r- j>d 2 y (toO = ~ jd 2 ye a D_<j = -i- Id 2 ye°£ia = ~£lA. (66) 
Therefore, if we choose the reference term P as 

P := ^£ T A, (67) 

where I is a future-directed out-going null vector field of a background reference spacetime, then P becomes, 

P(u, v) - ~£iA + -^-£jA. (68) 
sir 8tt 

Thus, the quasi-local integral P(u,v), which is to be interpreted as the quasi- local linear momentum, is given by the 
rate of change of the area of a given two-surface along the future-directed out-going null vector field I relative to I. 
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C. Quasi-local angular momentum 



Let us also write down the two-surface integral (|46|) in a geometrical way. Notice that the Lie bracket of the two 
null vector fields {n, 1} is given by 

[n, l] L = -F + «d a + {D_h)l. (69) 

Thus, ([46j) becomes 



l —j>d 2 ye a Un, (70) 
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where 

£a := ab e = e>afcC b . (71) 

This shows that @ is an invariant two-surface integral of the Lie-bracket [n, Z]l projected to the spacelike vector 
field £ :— £, a d a - If we choose the reference term L as 

L:=^-L 2 ye^ a [n,l]l, (72) 



then (Eq) becomes 



L(u,v,0 = ~ f<Pye°Mn, *]£ + ^ fl"<,<"U"- HI <" + i" <»=• ' ~ : '>*' 



It must be stressed £ a is an arbitrary function of {v,y b }. In particular, it need not satisfy any Killing's equations 
associated with isometries of a spacetime, or of a two-surface within a given spacetime. The quantity L(u, u;£), a 
linear functional of £ a , can be interpreted as the quasi-local angular momentum of a two-surface associated with an 
arbitrary function £ a (of y b ), as we will see later. 

D. Quasi-local Carter's constant 

Likewise, the fourth integral (|5^) can be written as 

J(u,v) = -±-j ) d 2 ye 2 °p ab A*[n, if + ±- jd 2 ye 2 ° p ab A*[n, l]l {d+A« = 0), (74) 

which may be interpreted as a quasi-local, finite, analog of the Carter's "fourth" constant, as we shall see in the next 
section. 

IV. ASYMPTOTICALLY FLAT LIMITS 

The equations ( |40| ) , (|4l]) , and (^2 ) turn out to be quasi- local energy, linear momentum, and angular momentum 
conservation equations, respectively! 1, and the equation ( |54"|) is interpreted as a quasi-local conservation equation of 
the generalized Carter's constant jl8, jl9L p0[ . In this section we shall evaluate the two-surface integrals (|65|), (|6"g|), 
((73|) and the associated flux integrals in the limiting asymptotically flat region where N2 — S2, and show that they all 
reduce to the well-known Bondi energy, linear momentum, angular momentum, and the corresponding flux integrals 
defined at the null infinity. The asymptotic form of the "fourth" integral ( [74] ) at the null infinity will be also computed, 
and it will be shown that it is proportional to the total angular momentum squared. 

In the limit where the afiine parameter v approaches to infinity, the asymptotic form of the Kerr metric becomes 

, , , , / 2m \ , 9 /4masin 2 i? 4rna 3 sin 2 ??cos 2 t9 \ , , 
ds 2 — > -2dudv - 1 1 )du z + - 1 \dudip 

9/ a 2 cos 2 i? \,„9 9 a 2 \,9 

+v 2 (l + —- 2 + • • • Jd-d 2 + v 2 sm 2 -&(l + — + ■■■ jdtp 2 

. 2n /4ma 3 8m 2 a 3 \ , , /a 2 sin 2 tf \ 2 

+ sm + —&r- + ■ ■ -) dvd v ~ {—jr- + ---)dv\ (75) 
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where d/du is asymptotic to the timelike Killing vector field at infinity. The asymptotic fall-off rates of the metric 
coefficients can be read off from the above metrical], ^2], ^3|, |24|, |25| ], 

e ff = w 2 (sini?){l + 0(4)}, (76) 
H)(l - 2(^+0(1,}, (78) 



Pipip 



P* v = 0{\), (79) 

2/1 = 1 +0(4), (80) 

Ajf = ^_ + 0(4 (81) 

^~ + 0(l), (82) 
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A ± ^0(-). (83) 



From these asymptotic behaviors, we can deduce the fall-off rates of the following derivatives, 

d+a = 0(\), d-a=-+0(\), d+p a b = 0{-), d-p ab = 0(\), £^p a b = 0(-), 



= -4wsim? + 0(1), ir a = ~2vsmti + 0(1), ir ab = --e a p ac p bd (d+p cd ) + 0(1), 



n v = 6masin 3 tf + O(-), tt^ = 0(4t). (84) 

A. The Bondi energy-loss relation 



Since the integrand of the r.h.s. of ( |49| ) assumes the typical form of energy-flux, we expect that it represents the 
energy-flux carried by gravitational radiation crossing S%. Then the l.h.s. of (pl9| ) should be the instantaneous rate 
of change in the gravitational energy of the region enclosed by S2. The energy-flux integral in general does not 
have a definite sign, since it includes the energy-flux carried by the in-coming as well as the out-going gravitational 
radiation. But in the asymptotically flat region, the energy-flux integral turns out to be negative-definite, representing 
the physical situation that there is no in-coming flux coming from the infinity. 

Let us now show that the equation ( |49| ) reduces to the Bondi energy-loss formula(l7| in the asymptotic region of 
asymptotically flat spacetimes. To show this, let us first calculate U(u,v) in the limit v — ► 00. Since the null vector 
fields n and I asymptotically approach to 

I — >cL, (85) 
the natural background spacetime is the flat spacetime so that the embedding degrees of freedom are given by 

A ± a = 0, 2h= 1. (86) 

That is, the background fields n and I become 

n = <9 + -i<9_, T=d— (87) 

Then it follows trivially that the total energy at the null infinity coincides with the Bondi energy Ub(u), 

lim U(u,v) := U B (u) = m, (88) 
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where m is the Bondi mass of asymptotically flat spacetimes. One can further show that the equation ( fl9| ) is just the 
Bondi energy-loss formula, 

d 



-U B (u) = - lim -J- (fdil v 2 p ab P cd (d + p ac )(d +Pbd ), (89) 
or, equivalcntly, 



4-u b (m) = i dn (d+cf < o, (90) 

du 167T Js 2 

where we used the expressions ([7?]) and (|7^) . Notice that the negative-definite energy-flux is a bilinear of the traceless 
current j a b defined as 

j a b := P ac d +Pbc (j a a = Q), (91) 
representing the shear degrees of freedom of gravitational radiation. 

B. The Bondi linear momentum and linear momentum-flux 



Let us now evaluate P(u, v) in ( |6q ) and the corresponding quasi-local momentum flux integral in the asymptotic 
region of asymptotically flat spacetimes. We find that the total linear momentum P(u, v) becomes zero in the 
asymptotic limit, 

lim P(u, v) := P B {u) = 0, (92) 

v— *oo 

from which wc infer that the total momentum flux is zero, 

^-Pb(u)=0. (93) 
du 

The result (^3|) can be also obtained by evaluating each term in the "Hamiltonian" function H (fll]). To evaluate 
the momentum- flux term by term, let us notice that the fourth and the seventh term in (|34|), which are non-zero 
individually, add up to zero asymptotically, 

^e-y abPcd n ac n bd + i-e" p ab p cd (D+p ac )(D +Pbd ) = p ab p ad (D_p ac ){D_p M ) 

where we used the definition (^9|) of ir ab . All other non- vanishing terms are given by 

lim J_ Id 2 y( -fur* id) = - (95) 
lim JL <fd 2 y (le-v-Kh^) = 1, (96) 
lim -i- (fd 2 ye< T R 2 = ] X , (97) 

where x = 2 for a two-sphere S 2 - Therefore we have 

4-Pb(u)=0. (98) 
du 

C. The Bondi angular momentum and angular momentum-flux 

The total angular momentum at the null infinity is naturally defined as the limiting value of the general quasi-local 
angular momentum L(u,v;£) in (|73|), 

lim L(u,»;0 :=L B (u;0. (99) 
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Since the background fields n and I in ( |87| ) commute, 

[n, I]l = 0, (100) 

it follows that 

£ B («;0 = o (ioi) 

for all £. Let £ be asymptotic to the azimuthal Killing vector field of the Kerr spacetime such that 

£:=£°d ^— • (102) 

Then, we have 

L B (u;£) = / dy? / di? (6ma) sin 3 z? 

16tt Jo 7o 

= ma, (103) 

which is just the total angular momentum of the Kerr spacetime. 

The total angular momentum flux at the instant u is given by the asymptotically limiting form of the the equation 
([42|) , which is 



= lim -L I d 2 y (7r ab £ iPab + 7r <T £ i a - h£^ h - A^£ i 7r a ). (104) 



Let us evaluate each term in the r.h.s. of this equation. The first term is given by 

7T ab £^a b = {-^e°p aC p bd (d +Pcd ) + 0(l)}£sPa b 



= -sini?(9 + C)(^C) + 0(i), (105) 

so that we have 

<f d 2 yn ab £(:p ab — ► - <f ffi {D+C){£ ( C). (106) 

The second term becomes 

it a £ i (j= I -2usim? + 0(l)}.%7. (107) 

Since cr is asymptotically given by 

cr = 21nv + In |siiM?| + ln|l + o(— )}, (108) 

we have 

so that the second term becomes 

j>d 2 yTT a £ i a = o(~) — ► 0. (110) 

The third term becomes 

hj£f7I7, = -7T h £(h + £((h7T h ) 

= £ i (4msm-d + hTr h ) (HI) 
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where we used that 



Thus we have 



The fourth term is of the order of 



so that 



£ 5 (sim9) = 0. 
j>d 2 yh£^ h =o(^j — >0. 

d 2 2/A + a £ 5 7r Q = 0(\) — »0. 



(112) 



(113) 



(114) 



(115) 



If we put together ( |l06| ) , ( |11C| ), (113), and flllq ) into (|104[) , then the total angular momentum flux at the null infinity 
is given by 



lim 



d 2 ye tT p ac p M {d + p cd ){£ i p ab ), 



or, equivalently, 



dL B 
du 



__L I dn {d + o(£ t c), 

!6tt Js 2 



(116) 



(117) 



which is precisely the Bondi angular momentum flux at the null infinity |2C 



D. Gravitational Carter's constant 



Let us define the asymptotic quantity Jb{u) as 



lim v 3 J(u,v) := Jb(u). 



Because of the equation (10C), the subtraction term Jb(u) is zero, 

J B (u)=0. 



Then the asymptotic integral Jb (u) becomes 



J B (u) = 2{maf 



(118) 



(119) 



(120) 



Thus, Jb(u) is (twice of) the angular momentum squared, deserveing the name the gravitational analog of the Carter's 
"fourth" constant at the null infinity. 



V. IN-GOING NULL COORDINATES 

One might be also interested in applying this formalism to black holes, and try to obtain quasi-local quantities 
defined on the black hole horizon and corresponding fluxes incident on that horizon. For this problem, it is appropriate 
to choose a coordinate system adapted to the m-going null geodesies. Such a coordinate system is described by the 
metric 



ds 2 = +2dudv - 2hdu 2 + e a pab (dy a + A^du + A^dv) (dy b + A b du + A b dv) 



(121) 
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In this coordinate system, the future-directed out-going and in-going null vector fields n' and I' are given by 

n':=d+ + hd-, (122) 
i' :=-£_, (123) 



respectively, which are normalized so that 



The inverse relation is given by 



< n', I' >= -1. (124) 



d+ = n' + hi', (125) 
d- = -I'. (126) 



If we repeat the same analysis as in the previous chapters using the metric (121), we obtain a new (but equivalent, of 
course) set of quasi-local conservation equations, which we present without derivations. The conjugate momenta are 
found to be 

7T h = -2e a D_a, (127) 
n a = -2c a D-h - 2hc a D-a -c a D+a, (128) 
7r a = e 2 °p ab F + b , (129) 

n ab = he°p ac p bd D„p cd + ±e°p ac p bd D +Pcd , (130) 
and the "Hamiltonian" function H is given by 



TJ — rrr- I U ^ — & ^ 2 — 2(7 „ „ I ^ _ — <T „ n „0,C „bd 

ti = --e 7r h 7r CT + -fte 717, - -e p a b P ab P cd7T n 



8_ 

du 

d i 

—P(u,v) = -—(l)d 2 yH, (133) 

(134) 



--Tr^+a - ^ ab D +Pab + ^. c °p- b p- d (D + p ac ){D +Pbd ) + c°R 2 . (131) 
A new set of quasi-local conservations equations are found to be 

-U(u, v) = ^ j>d 2 y (Tr ab d +Pab + ir a d+<j - hd+n h - A«d+n a ) , (132) 
1 

n -P(u,v) = 

OU 167T 

^-L{u,v-0 = -^-jd 2 y (Tr ab £ (Pab + Tr a £ s cT - h£^ h - A + a £^ a ) {d+C = 0), 

= ^ fd 2 y(ir ab £ A+ p ab + 7r a £ A+ a-h£ A+ ir h ) (d+A« = Q), (135) 
where U(u,v), P(u,v), L(u, v;£), and J{u,v) are defined as 

U(u, v) := ^ I d 2 y (h tt h - 2e ff D+a) + U, (136) 
P(u,v):=^Id 2 y(ir h )+P, (137) 

L{u,v^):=^jd 2 y{e^ a )+h (138) 
J(u, v) := ^"d 2 2/ (A + Q 7r a ) + J, (139) 

wher e [/, P, £, and J are undetermined reference terms as before. Notice that we could also have written the equation 
(IH) as 

^f{u, v) = ^- j>d 2 y (w ab D +Pab + n a D + a - hD + -K^j , (140) 
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as in (|40|). In geometrical terms, these quasi- local quantities can be expressed as, 
U(u,v) := - — £ n/ A + —£ W A, 

07T OTT 

P{u,v) = ^-£ v A-^-£ v A, 

OTT OTT 

L(u,v;0 = ±-jd 2 ye°Un, l]l - ^ <f d 2 y e^ a [n' , I'M (d+C = 0), 

J(u,v) = ±- fd 2 ye 2 °p ab A«[n, l t - ^ f d 2 y e 2 ° p ab A + a [n' , Y]l (d + A« = 0). (141) 

Here n' , I' are future-directed in-going and out-going null vector fields of a background reference spacetimc 

ds 2 = +2dudv - 2hdu 2 + c a p ab (dy a + A^du + A^dv) (dy b + A*du + A}dv) , (142) 

such that 

d t „ d \ . r / d T „ d 



i 1 := ( — - A a — \ + hi— - A_ a — ) 
\du + dy a J \dv dy a I 



dy a / \dv dy a J ' 

dv dy a , 

VI. QUASI-LOCAL HORIZON 

Recall that the event horizon is a global concept which is inseparable from the notion of infinity. Therefore, in 
order to discuss the dynamics of black holes quasi-locally, we have to introduce a new notion of quasi-local horizon, 
which refers to the quasi-local region only. We define the quasi-local horizon H as a three-dimensional hypersurface 
on which the vector field <9+, which is an arbitrary vector field except that it is asymptotic to the timelike Killing 
vector at the infinity, has a zero norm, 

< d + ,d+ >= -2h = 0. (144) 
The location of the quasi-local horizon TL can be found by solving the equation 

h(u,v,y a ) = (145) 
for v, and the generator of the quasi-local horizon is given by 

which is out-going null on 7i. 

Let us remark a few properties of this quasi- local horizon. First, the quasi-local horizon is defined for generic 
spacetimes that do not have isometries in general, and its location is not fixed, but varies as much as the choice of the 
vector field 9+ does. In this sense the quasi-local horizon is not a spacetime invariant but a covariant notion. Even 
the signature of the vector field 9+ is not determined a priori. In the region where 9+ is non-spacelike, however, the 
quasi-local horizon may be regarded as a generalization of the Killing horizon, since it is defined as the hypersurface 
where 9+ has a zero norm. For instance, for the Schwarzschild solution, we have 

2h =l-—, (147) 

v 

so that h — on the Killing horizon v — 2m. Moreover, since the quasi-local horizon is generated by the out-going 
null vector fields, all the fluxes crossing the quasi-local horizon is purely in-going, just like the stretched horizon of 
Price and Thorne. 

In this section, we shall delimit our discussions of the quasi-local conservation equations to the quasi-local horizon 
7i, and find that the quasi-local conservation equations on TL coincides exactly with the quasi-local conservation 
equations of Price and Thorne pffl defined on the stretched horizon. 
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A. Surface gravity « 

In order to discuss the dynamics of quasi-local horizon, it is useful to introduce the notion of the surface gravity k 
to the generic, time-dependent, quasi-local horizon. On the quasi-local horizon TL on which h = 0, the vector field x 
defined as 



X--=-- A +l^. (148) 



du + dy a 
becomes a generator of TL, since we have 

X • X = -2h\ n = 0. (149) 
Hence Va(x ' x)\h is normal to TL, which means that there exists a function k defined on TL such that 

Va(x-X)\h--=-2kXa\h, (150) 

where A is a spacetime index such that A = {+,— ,a}. Notice that this function k can be defined on any null 
hypersurface. When the null hypersurface coincides with the event horizon, this function is the surface gravity of 
the corresponding black hole. But one may use the same terminology for k on a quasi-local horizon TL, since it is a 
(segment of) null hypersurface beyond which a local observer whose worldline 9+ has a zero norm on TL does not have 
access to. It may be instructive to notice that, for the Kerr black hole, x is the generator of the event horizon, since 
it becomes 

x-Tu-^h < 151 > 

where d/du and d/dip are timelike and axial Killing vector fields, and fin is the angular velocity of the Kerr horizon 
relative to the Killing time. 

Let us compute n on the quasi-local horiozn TL. In the basis {d/du, d/dv, d/dy a }, the components of x are given 

by 

X A = (1, 0, -A«), 

XA := 9ABX B = (-2fc, 1, 0). (152) 



If we put (152) into (15C), then we find that 

d+h\ H = d a h\ H = 0, (153) 

and that K is given by 

k = D-h\ H . (154) 
Notice that, in general, k is not constant over TL so that 

d+K ^ 0, d a K ^ 0, (155) 
which reflects the dynamical nature of the quasi-local horizon TL. 



B. Quasi-local energy conservation on TL 



Notice that if we restrict the energy equation (140) to the quasi-local horizon TL, then it becomes, 
du i on j n 



-Un = ^i d 2 y \le°p ab p cd (D +Pac )(D + p bd ) - e°(D + a) 2 - 2e^D + a}, 
i lbn J n 1 2 J 



1 



U n := -— <b d z ye a D + a + U H . 
This equation is identical to the integral of the following equation 



(156) 
(157) 



du 



8tt 



16tt 



H 



8tt 



M ab 
J ab a H 



(158) 
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over the stretched horizon of Price and Thorne, where the notations are such that 

6 n = D+a, 

w 1 * 

<b = 2 C ^+^6, 

^.afc A.acihd'H ^ n — a ac bd t-i „ 

<r w = 4> 4> o cd = -e p p V + p cd , 

K = D-h\w (159) 

This equation was studied in detail in Eq. (6.112,E) in |27j ]. 

It is interesting to discuss the limiting case when the quasi-local horizon TC coincides with the event horizon. When 
this happens, the area A-h of TL always increases due to the area theorem, so that we have 

^ = / d 2 y (e a D+a) > 0. (160) 
du J H 



Furthermore, if the subtraction term U-h is chosen zero, then by the equation (157), U-h is non-positive, and when the 
black hole no longer expands so that D + a\^ = 0, then U-h becomes zero. For instance, for a Schwarzschild or Kerr 
black hole[^7|, we have 

U H = if U H = 0. (161) 

This counter-intuitive aspect is a manifestation of the well-known teleological nature of the event horizon. That is, 
when the event horizon evolves, its quasi-local energy must be negative so as to cancel out the positive in-flux of energy 
carried by subsequently in-falling matter or gravitational radiation, leaving Uu — when the black hole reaches the 
final stationary state. 



C. Quasi-local momentum conservation equation on "H 



Let us mention that the momentum equation (133) has a similar structure to the integrated Navier-Stokes equation 
for a viscous fluid |E8|, 



dPi, 
du 



dS k (pSik + pviV k - a' ik ^j , 



where P; and a' ik are the total momentum and the viscous term, 



(162) 



E t = dV{ P Vi), 



'ik 



1) 



^dx k 



dvk 



2 dm\ 



(Sik 



dvi 
dx l ' 



(163) 
(164) 



and 77 and C are the coefficients of shear and bulk viscosity, respectively. This equation tells us that the rate of the net 
momentum change of a fluid within a given volume is determined b y th e net momentum-flux across the two-surface 
enclosing the volume. Notice that the "Hamiltonian" function H in (131), which is at most quadratic in the conjugate 
momenta 717, assumes the form of momentum-flux of a viscous fluid. Namely, terms quadratic in 717 may be viewed 
as responsible for direct momentum transfer, terms linear in 717 as viscosity terms, and terms independent of tti 
as pressure terms. From this point of view, one may interpret the "Hamiltonian" function H as the gravitational 
momentum-flux and the two-surface integral 



Ph = 



1 

167T 



d 2 y{7T h ) + P 



(165) 



as the quasi-local gravitational momentum within N2. On the quasi-local horizon 7i, the equation (133) becomes, 



-2(7 ab 



P n a ir b - -e ir a ir h - ~-K h D + o 



(166) 
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The first term on the r.h.s. is given by the Euler number 

1 



16tt 



d 2 ye a R 2 



7 X) 



(167) 



H 



where x = 2 for a two-sphere. The second and third terms are quadratic in the momenta, and the last term is linear 
in the momentum. It is curious that this conservation equation is missing in the work of Price and Thornc p7| . Let 
us note that, in terms of the configuration variables, the integrand of the r.h.s. of (166) can be written as 



(168) 



D. Quasi- local angular momentum conservation equation on Ti 



In this section, we shall show that the equation ( 134 ) when restricted to the surface Ti coincides with the quasi- loca l 
angular momentum equation of Price and Thorne on the stretched horizon. Let us first notice that the equation (134) 
on Ti can be written as 



0_ 

du 



Ln = lo7r { ~ C,T(2K + D + a ) £ i a + \^p ac p hd {D +Pab ){£ iPcd ) - ^ + a ^7r Q }, (169) 

(170) 



L ^ : =T7Tf d 2 y(C a n a ) + L H 



The equation (169) turns out to be identical to the angular momentum conservation equation of Price and Thorne, 
which is given by 



H 



1 



1 

loTT 



n 



H 



16tt 



e^PabF, 



(171) 
(172) 



To show this, let us write the equations ( |L69| ) and (17C) as 

lL Ln = lb f H d2y { 2e<7£ Z K + e °M D +<*) + \c a p ab £d D +Pab) + e 2 ° Pab F + «£sA^, (173) 

Lu := lb £ ( e2a P" bF +-^ + Z «' ( 174 ) 
where we used the identity 

/ d 2 y£t:f= I d 2 yd a (Cf)=Q (175) 
Jn Jn 



for any scalar density / with the weight — 1. Using the definitions of Lu in (174) and 11^ in (172), we obtain the 
following identity, 

= - (f d 2 y{e°e(D + a)nK + e°(D + nn a H +e°eD + TlZ + e°TlZ{A + , t] a A 
= - ^d 2 y{e^ a D + n o H +^e(D + a)U^+c^( K D + e + [A + , flg)} 
= -<f d 2 ye°e(D + I% + d H I%), 



(176) 



where we used the diffiV2-covariant derivative of £°, 

D + e := o+e - [a+, e]£ 



(177) 
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defined in the section |, and the condition that 

d+C = 0. (178) 
The first and second term on the r.h.s. of ( |173p become, 

±- f *y + *£<{D + *)} = f/ye^H^ +^°n, a ) , (179) 



and the third term on the r.h.s. of (173) is given by 

±e°p ab £ e (D +Pab ) = \e a p ah {eV c (D +Pab ) + (D +Pc6 )(V„0 + (D + p ac )(V b {; c ) - (V c ^)(D + p ab )} 

= \o^ c {v c {p ab D +Pab ) - {V c p ab )D+p ab ) + c°p ab (D+p bc )(V a O 

= e°p ab {D +PbL ){V a O 

= -Cc°p bc V b (D + p ac ) + e°p ab V a (eD +Pbc ) 

= -2Cc°a Ha b ; b + V a (ee a p ab D +Pbc ). (180) 
Here we used the unimodular condition (|3l]), and the metricity condition 

V a a - V aPbc = 0. (181) 



Notice that the shear tensor a Ha is given by 



i n b := <p bc a nac = \p bc D +Pac . (182) 



Therefore we have 



/ d 2 yh°p ab £t(D + p ab ) =-<f d 2 y2Ce°a Ha b ]b . 



(183) 



If we put together (f79), (179), and (183), then (173) becomes 



^d 2 ye°c{D + I%+8 n nK + ±- K , a+ -L-d n ,a - ^na fi} = 0, (184) 



for an arbitrary function £° that satisfies (178). This shows that the two equations (|169|) and (|17l|) are identical 



VII. SUMMARY 

In this note, I presented a set of quasi-local conservation equations that were found while studying the Einstein's 
equations using the (l+l)-dimcnsional description, and studied physical significances of these conservation equations. 
The key observation is that the affine coordinate v can be treated as a natural time coordinate in this formalism. The 
subsequent Hamiltonian formalism was developed with respect to this time coordinate. One of the outcomes of this 
analysis is the discovery of a set of Bondi-like quasi-local conservation equations for the vacuum general relativity, 
which reproduce both the well-known conservation equations in the asymptotic null infinity in asymptotically flat 
spacetimes and the corresponding conservation equations on the inner quasi-local horizon of a generic dynamical 
spacetime. All of these quasi-local quantities are expressed in geometrically invariant terms such as the area of the 
two-surface and a pair of null vector fields orthogonal to that surface. It was also found that each quantity has 
a natural interpretation as the quasi-local energy, linear momentum, and angular momentum of a two-surface and 
corresponding fluxes crossing that surface. 

In addition to the above quasi-local quantities, we also obtained the quasi-local analog of the Carter's "fourth" 
constant of gravitational field, which is somewhat like the angular momentum squared, measuring the "intrinsic" 
angular momentum of a two-surface. The Carter's "fourth" constant is known to exist for spacetimes that possess 
two commmuting Killing vector fields such as the Kerr black hole. But in our analysis, it was found that the quasi-local 
analog of Carter's constant exists under the condition 

d+A^ = 0, (185) 
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which is much less restrictive than the existence of two-commuting Killing fields. 

Dynamics of the quasi- local horizon was discussed briefly, but it deserves further studies. Applications of this 
formalism to astrophysical problems involving black holes and gravitational radiations are extremely challenging. 
These problems are left for future works. 
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